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I. INTRODUCTION 



The pion form factor, Ft^{s), in the p resonance region has been the subject of several 
recent experiments [l|-l6[ which, together, provide the most precise picture yet of the 
structure of the p resonance, including the prominent p — u interference 'shoulder'. The 
markedly improved precision of the new data (e.g. a statistical uncertainty of 1% versus 
4-5% for earlier studies j3"10|) enables high-precision studies of isospin breaking in the 
p — uj sector and an improved determination of the dominant tttt contribution to the 
leading order hadronic vacuum polarization term in Standard Model (SM) estimates of 
the anomalous magnetic moment of the muon, a^ = {g — 2)^/2. 

Isospin breaking in the p — u sector is of particular and ong oing importance to studies 
of charge symmetry breaking in the A^A^ interaction (see Ref. and references therein 
for more details). Significant unresolved questions include i) the precise scale of the 
p — UJ mixing matrix element, Ilpu) ~ {p\H\uj), in the region of the p and ii) the scale 
of the contribution associated with the unavoidable direct (i.e. non-mixing induced) 
isospin-breaking w — )■ vrvr transition, which has often been neglected in determining IIp^. 



Previous analyses of these issues (e.g. in Ref. 12|) were hindered by data of relatively 
lower precision, making a precise determination of the phase between the p and u con- 
tributions impossible. Without a well- constrained phase, it is not possible to obtain 
a well-constrained estimate of the direct u — >■ tttt contribution, and the uncertainty in 
IIp^ (m^) remains large (in fact much greater than what is often quoted, for example in 
Ref. fllj). The recent wealth of high-precision e"'"e~ data in the resonance region repre- 
sents a good first step in addressing these concerns. In this paper we revisit p — uj mixing 
in Ft^{s) and use the recent data to effect a separation of w — )■ vrvr and mixing contribu- 
tions, from which we determine the ratio of isospin-breaking (IB) to isospin-conserving 
(IC) couplings, gui-KTrl Qpi-KTT) and the off-diagonal vector meson self-energy matrix element 
np^(m^). We perform this analysis using a selection of models for Ft^{s) in order to clarify 
the numerical significance of the unavoidable model dependence of such a separation. 

One important reason for the recent interest in F-,j{s) is the dominant role it plays 
in theoretical estimates of the leading order (LO) hadronic contribution to a^, [a^]^^- 

is now known experimentally to a remarkable 0.5 parts-per-million precision as a 
result of BNL experiment E821 isf. It is well known that [a^had which, after the purely 
leptonic contribution, is the largest of the SM contributions to a^, can be obtained 
from a weighted integral of the e^e~ — ?• hadrons cross-section with the dominant 
contribution coming from the broad p resonance. Precise data on the e^e~ — )■ vr+vr" 
cross-section (or, equivalently, -^^(s)) below 1 GeV is, therefore, a crucial ingredient in 
the determination of [a,p\had- Use of recent e+e^ — t- tt+tt^ data as input for the tttt part of 
[a^had leads to SM estimates for which are consistent with one another (even though 
the form factor measurements themselves do not all agree on the precise s-dependence 
of -F^(s) [lH), but which deviate from the experimental value by ~ 3.5 - 3.8(T 



(It should be noted, however, that preliminary radiative return results from BaBar differ 
from previous electroproduction results (22|.) 

As is well known j23-27|, in the isospin limit, CVC allows the s < isovector 
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contribution to [o//]^^ to be computed using hadronic r decay data 28l-[33| in place of the 
corresponding isovector electroproduction (EM) data. At the present level of precision, 
additional small isospin-breaking corrections to the CVC relation 21, 23-261 33"40| must 
be taken explicitly into account. We denote the collection of such corrections by [f^a^]^^ 
in what follows. Incorporating the isovector r decay data in this manner, one obtains 
estimates for the SM contribution to which deviate from the experiment by ~ 1 — 
2(7 HHHiil. The low-energy region from which the bulk of the SM contribution to 
\a^]j^2i is generated is dominated by the vrvr component of the EM cross-section. The IB 

are 



21, 36-40 



corrections in this region, which have been extensively studied in Refs. 
believed to be well understood. However, the component of these corrections associated 
with p — LO mixing, which in general is model dependent, has, to date, usually been 
estimated from fits to Ft,{s) which employ only a single model for F^(s) (usually that 
of Refs. 36, 13]). This procedure turns out to lead to a significant underestimate of the 
resulting uncertainty, and hence also of the uncertainty in the r decay-based value of the 
two-pion contribution to 



^LO 

''t^ihad- 



In Ref. j4l| several models were used to examine the model-dependence of [6a 



f^ihad- 



It 

was found that the variation across models was significantly greater than the uncertainty 
(associated with errors in the experimental data) obtained for any individual model, 
though not to such an extent that it would resolve the r vs. e~^e~ discrepancy. It should 
be borne in mind in this regard that finite energy sum rule studies of the electromagnetic 
current-current correlator j42|, based on the then-current EM database, yielded values of 
as{Mz) ~ 0.1140 — 0.1150, which are ~ 1.5 — 2a low compared to recent high-precision 
lattice determinations 43|. In contrast, analogous studies of the charged isovector current 
correlator ji^] , based on the final version of the ALEPH hadronic r decay data j3l[ , yield 
a value of as{Mz), in excellent agreement with the new lattice results. The preliminary 
BaBar results for the e~^e~ — )■ tt+tt" cross-sections reported by Davier at Tau'08 |22l] . 
moreover, are in better agreement with r expectations than with previous EM results, 
also favoring the r determination, though the r-EM discrepancy has been reduced by the 



recent re-assessment of IB corrections performed in Ref. [2l| . The question of whether the 



determination of [a^]^^ incorporating r data, or that based solely on EM cross-sections, 
is the most reliable thus remains an open one. The goal of achieving as accurate as 
possible an understanding of the uncertainty on the IB correction, [Sa^]f^^^, to be applied 
to the raw r-based version of [o^]^^ thus remains a highly relevant one. 

In the present paper, we update the analysis of the p—u mixing contribution to [5a^]^^ 
contained in Ref. 
data 6! 



41 



—7- TT^TT' 



taking into account the corrected version of the SND e+e" 
the new CMD-2 e^e~ — )■ tt+tt" results j2|, and the latest Initial State Radiation 
(ISR) determination of F^^ls), from KLOE In addition, we provide details of the fit 
results associated with the various models of -^^(s), and new results for the separation 
of the IB part of the EM cross-sections into pure p — u mixing and direct w — )■ tttt 
coupling induced contributions. We find that, taking the model dependence of the p — u 
contribution to [^O/i]^^ into account, the uncertainty assigned to this quantity should be 
increased by 1.5 x 10"^''. The improved data also allows for a more reliable separation 
oi p — u mixing and direct non-mixing — > vrvr contributions to Ft^{s). The uncertainty 
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in the scale of the non-mixing contribution, while still large due to weak constraints on 
the relative phase of the uj and p contributions to the amplitude, is reduced by a factor 
of two compared to the analysis of Ref. [l2| and favors a non-zero value. As a result, 
the central value of Iipui{m?p) is significantly different from that obtained using analyses 
in which direct uj contributions are ignored. 

The rest of this paper is organized as follows: Section [TTl reviews the different models 
of -F^(s) considered in our analysis; Section UTTl discusses the four recent e"'"e~ — )■ vr+vr" 
data sets and the results of fits to each set using the models of Section [TTl Section [IV] 
updates estimates of the p — u contribution to [^a^]^^; and, finally. Section IVl reviews the 
(model-dependent) separation of the direct a; — t- tttt and mixing contributions to Ft^{s) 
and presents the results of that analysis for the models and data sets of Sections [TTl and 

m 



II. MODELS OF F^{s) IN THE RESONANCE REGION 

Many recent estimates of the p — u mixing contribution to [f^o^]^^, which we denote 
by [Scitj]had-mix what follows, have often been based on a model of Ft,{s) which supple- 



ments the ChPT-based formulation of Guerrero and Pich (GP) [45[ with a p — w mixing 



contribution as proposed in Refs. |36l. l37l| (GEN). We refer to this model as the GP/GEN 
model. However, -F^(s) is well-described by many other models in the literature. Various 
analyses have used the Kiihn-Santamaria (KS) model ji^, the Hidden Local Symme- 
try (HLS) model [13], and the Gounaris-Sakurai (GS) model [il]. Each model differs 
in the manner in which the dominant p contribution to -Ft(s) is represented, with the 
GS and GP/GEN models explicitly taking account of strong final state interactions in 
e'^e'^ — )■ TT+TT^. The use of models for -^^(s) is unavoidable at present if one wishes to 
separate the IG and IB contributions to the EM cross-sections in the interference region, 
this separation being necessary for an estimate of [ScL^l]had■m^x■ 

The KS model [i^ combines contributions from the p, p', and p" resonances in such a 
way as to satisfy the charge constraint F^(s = 0) = 1. An u contribution of the nominal 
p — u mixing form is added to account for the narrow interference shoulder. The resulting 
expression for is 4^ 



^(KS)(^ 



1 + /3 + 7 

where the vector meson propagator is given by 



my — s — irnvTvi-s, my, T\ 



Pvis) = ... . JT .... . . (2) 



and 6, (3 and 7 are complex-valued constants. Unless otherwise stated, Tv{s,mv,Tv) 
is the standard s-dependent width implied by p-wave phase space for the decay of the 
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vector meson V, explicitly for the p, 



Tp{s,mp,Tp) 



m. 



(3) 



where (3{s) 



1 _ and 

s ' 



PJs) 



m. 



(4) 



{ml - s - im^r^) 

with approximated by a constant. 

An alternate version of the KS model, referred to in what follows as KS', is also 
sometimes employed in the literature. The KS' form for -F^(s) is 



^(KS')( 



'Pp{s) (l + 5-|P^(s)) +/3P,,(s)+7Pp"( 
1 + /3 + 7 



(5) 



It is this alternate version which is employed in Ref. j2l|]. Because of strong cancella- 
tions in the region of the p — u interference shoulder, the KS and KS' models produce 
significantly different results for [^CLfilhad mix- return to this point below. 
The HLS model 47, 49 of -F^fs) has the form 



(6) 



with ttHLS a constant associated with a non-resonant 77r7r coupling. Despite having no 
explicit contribution from higher resonances such as the p', the model provides a good 
quality fit to the data below 1 GeV. It also turns out to produce the correct -F^(s) phases 
in the elastic tttt scattering region after the model parameters have been fitted [49i] . 
though the phase constraint is not imposed in the basic structure of the model. 
The GS model 48|, as adopted by CMD-2 [1] and used here, is given by 



BWf^)(s) 



where 



BW 



(GS). 



l + 6—P^is) 



(l + d{my)^ 



(7) 



[m 



with 



d{mv) 
H{s) 



3m^ 



V 



in 



s + f{s, my, Fy) — imvTv{s, my, Fy)) 



(my + 2j9^(m^)) 



2m^ 



+ 



my 



m^mv 



dH 



pi{s)[H{s) - H{m'y)] + {mi, - s)pi{mi) — {mi\ 



(9) 
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where Ptt{s) = — is the pion CM momentum for squared invariant mass s and 
Ty = Ty{my, niy, Ty). The elastic vrvr scattering phase constraint on F^(s) is reproduced 
by the GS model by explicit construction. 

For the KS, HLS, and GS models, the constant 6, which parameterizes the strength 
of the narrow IB amplitude, is taken to be complex since a non-zero phase is, in general, 
unavoidable in the presence of an IB direct (non- mixing) — t- vrvr contribution |12| . 

Resonance chiral effective theory provides the basis for the GP model of Ft^{s), which 
is given by |45|] 



Ff^)(.) = P,(.) exp 



—s 



967r2/2 



ReL( ^ 



/ 2 2 

/ ml mt 



mt 



+ ^ReL 



' 2 

s mi 



(10) 



where 



2 2 



L\—,^\=in\^\ + 



m 



5m 



S 6 



:iii 



"p/ \ ■ ■ "p / 

with (5{s) = a/1 — Am^ / s and the s-dependent width, T p{s,mp,T p) appearing in Pp{s) 
(Eq. [2]) replaced by the resonance chiral effective theory expression 



rp(3) 



mpS 



- (eis - Ami) PAs)' + ^Ois - Am\) /3k(s)') • 



:i2) 



Because the GP model implements the constraints of chiral symmetry explicitly at next- 
to- leading order in the chiral expansion, it also ensures the correct phase for F-j^i^s) in 
the elastic region to this same order. In Ref. 3 

appearing on the RHS of Eq. f[T^ is allowed in order to account for the 



a small rescaling of the coefficient 

1.5 MeV 



contribution of 7r7r7 decays to the total width of the p [37 
rescaling by multiplying the RHS of Eq. ([T^ by 



50[. We implement such a 



1 + 



rp(?n2) 



(13) 



with bV p a fit parameter. Some IB effects are incorporated into Fn'~'^\s) if one evaluates 
the phase space factors in the s-dependent width using the physical charged vr and K 
masses. The GEN modification of Fi'^^\s), designed to incorporate the p — u mixing 
contribution not included in the original GP model, has the form 



FfP/OE^)(s) = Ff''\s) 



Ppis) 



e 



3m^ 



-I PUs) . 



(14) 



The parameter 9p^ 
As discussed in Ref. 



37 



was assumed real in Ref. 

|4l| , the original version of the GP / GEN model does not provide 
a good fit to the corrected version of the earlier GMD-2 data [ij. As already noted, the 
fact that the p — u "mixing" signal is actually a combination of mixing and direct a; — )■ tttt 
effects means that an effective representation of this combination using the form given 
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by Eq. ( 1141) is not generally possible without allowing 9p^^ to have a non-zero phase |12[ . 
In Ref. [41] the GP/CEN model was extended in this way, treating the phase of 6p^ as a 
fourth parameter to be fit to the data, and an acceptable fit became possible. We refer 
to this version of the GP/CEN model as GP/CEN"*" below. In most of what follows, 
however, we will focus instead on an alternate version of the GP/CEN model, referred 
to below as GP/CEN^"*", which incorporates the effects of strong tttt rescattering not 
only in the isospin-conserving e"*"e~ — > p — )■ 7r"'"7r~ contribution (the first term on the 
RHS of Eq. m]), but also in the p propagator factor of the isospin-breaking p — uj mixing 
contribution, leading to 



,s) . (15) 



where, as in the CP / CEN"*" version, Op^i is allowed to have a non-zero phase. 



III. INPUT DATA AND FITS 

Recent information on Ft^{s) in the interference region has been provided by the CMD- 
2 and SND experiments at Novosibirsk and KLOE at DAFNE. Results are available from 
a total of five separate data-taking runs. Preliminary ISR results from BaBar have also 
been reported, but are not yet publicly available. 

CMD-2 has reported data on e^e~ — )■ vr"*"-/?" in the centre of mass energy range from 
0.6 to 1.0 GeV, taken in two separate runs, one in 1994-95 [l| and one in 1998 [2]. We 
refer to these as CMD-2 (94) and CMD-2 (98), respectively. The 1998 run has five times 
the luminosity of the earlier run, and correspondingly smaller statistical errors, but a less 
precise beam energy calibration. Whereas the beam energy was measured very precisely 
using resonant depolarization for CMD-2 (94), this technique was not available for CMD- 
2(98) |5l| with the result that the beam energy for the latter is only known to a few parts 
in 10~^ [2] ■ A significant part of the beam energy uncertainty, moreover, is fully correlated 
(5l| . The combination of increased statistics but poorer beam energy calibration resulted 
in no overall improvement in the determination of the tttt contribution to [di^had from 
CMD-2(98) as compared to CMD-2(94)[2]. Indeed, the beam energy contribution to the 
systematic error in [a^]^f^ increased from 0.1% in CMD-2(94) to 0.3% in CMD-2(98)0. 
The less precise energy calibration of CMD-2 (98) also has an impact on studies of isospin 
breaking in -F-^(s) because such studies require a precise experimental characterization of 
the interference shoulder. 

SND results for the e^e' — t- tt+tt^ cross-sections are based on data taken from 1996- 
2000 with centre of mass energies in the range 0.39 to 0.97 GeV [sl, 0]. The SND data 
carries a slightly larger systematic error than CMD-2 (94), but much smaller statistical 
uncertainty and slightly better beam energy calibration. The CMD-2 and SND data sets 
are in good overall agreement. 

The most recent, KLOE, measurements of the e^e~ — )■ tt+tt" cross-section are ob- 
tained using the ISR technique [sl, 0]. It is well known that the earlier, 2004, KLOE 
results for Ft,{s) based on data collected in 2001, deviated in shape from both the 
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CMD-2 and SND results, lying systematically higher than either below ^/s = 700 MeV, 
and systematically lower beyond this [l, [sj. We refer to these results as KLOE(Ol) be- 
low. A new analysis has recently been completed using data collected in 2002 and an 
improved analysis strategy The new data represents twice as many events as before 
and corresponds to an integrated luminosity that is 1.7 times greater, and the analysis 
yields significantly reduced statistical and systematic uncertainties, with the systematic 
uncertainties now dominant. In addition, part of the previous shape discrepancy has 
been resolved, KLOE now being in better agreement with CMD-2 (98) below the p peak. 
A reduced disagreement persists above the p peak. In view of the fact that the results 
of the new analysis supercede those of KLOE(Ol) our study of IB in the interference 
region, reported below, is based on the former, which we refer to for clarity as KLOE(02) 
in what follows. 



Our earlier analysis of [^o,fi]had-Tnix^ reported in Ref. 41|, was based on the CMD-2(94) 



and uncorrected SND data sets. The KLOE(Ol) results were also used, though, without 
full covariance information, the resulting fits were unreliable. The fits presented below 
reflect the new CMD-2 (98) data, the corrected SND data, and the new KLOE(02) data, 
including full covariance information, where available. 

We have performed fits to each data set using the models described previously in 
Section |Tll In all cases the results shown below correspond to the bare form factor. 



that is the form factor with the effects of vacuum polarization removed |52l . |53| . Where 



applicable, the following input values have been used: m^; = 782.68 MeV, T^^ = 8.68 



MeV, rup' = 1465 MeV [55|, Tp/ = 400 MeV (310 MeV for the GS model ^]), A = 92.4 
MeV, = 139.57 MeV, and mx = 495 MeV. Although a contribution from the p" 
appears in the KS model, we have found that, because the data we consider is limited 
to ^/s below 1 GeV, it provides insufficient constraints on this term. We manually set 
the p" contribution to zero in the KS model to avoid it being used by the minimization 
algorithm to 'fine tune' the p region of the fit without regard for the consequences at 
higher ^/s. For similar reasons the coefficient of the p' contribution, /3, is taken to be 
real. Fit parameters and their covariances were obtained using MINUIT. The results for 
each model and for each data set are shown in Tables [U [Tll llllj and IIVI A blank entry 
indicates that a fit parameter is inapplicable to that particular model. For the GP/CEN"*" 
and GP/CEN+"'' models, the effective value of Tp is shown in brackets to highlight that 
it is in fact STp which is the fit parameter. Note also that for both of these models the 
quantity 16*^^1 is replaced by \6\ = |6'p^|/3m^ in the tables to facilitate comparison with 
the other models. 

It is apparent from Tables [H [III HHj and [IV] that, for each given model, all four data 
sets are in broad agreement with each other. The different models are also in general 
agreement with each other for each given data set. One exception is the GP/CEN"*" model, 
which yields a value of |5| ~ 20% larger than in the other models. This deviation, however, 
disappears when the effects of rescattering are incorporated into the p propagator factor 
which enters the p — u interference contribution through use of the modified GP/CEN++ 
form. Table [IV] reveals that fits to the KLOE data stray most from the fits to the 
other three data sets, in particular through substantially different central values for the 
magnitude and phase of the isospin-breaking parameter 6, and a much greater model- 
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TABLE I: Results of fits to the CMD-2(94) data. 



Parameter 


KS 


HLS 


GS 


GP/CEN+ 


GP/CEN++ 


rup (MeV) 


772.42±0.63 


773.84±0.62 


774.67±0.64 


775.89±0.61 


775.89±0.61 


Tp (MeV) 


140.47±1.34 


143.89±1.52 


144.37± 1.43 


(145.51) 


(145.47) 


6Tp (MeV) 








-1.83±0.63 


-1.86± 0.63 


\6\ (10-3) 


1.77±0.14 


1.84±0.14 


1.84±0.14 


2.25±0.17 


1.91± 0.14 


Arg((5) (deg) 


14.9 ± 3.3 


14.1 ± 3.2 


14.3 ± 3.2 


14.1 ± 2.9 


14.2 ± 2.9 


/5 


-0.136±0.005 




-0.073±0.005 




- 


(IHLS 




2.3479±0.016 






- 


xVdof 


37.1/38 


36.6/38 


35.8/38 


40.5/39 


40.6/39 




TABLE IL Results of fits 


to the CMD-2(98) data. 




Parameter 


KS 


HLS 


GS 


GP/CEN+ 


GP/CEN++ 


rup (MeV) 


772.65±0.55 


774.07±0.54 


774.87±0.55 


776.86±0.51 


776.86 ± 0.52 


Tp (MeV) 


141.56±0.81 


146.51±1.01 


146.61± 0.91 


(147.66) 


(147.63) 


(JPp (MeV) 








-0.32±0.56 


-0.34± 0.56 


\5\ (10-3) 


1.88±0.068 


1.94±0.069 


1.93±0.069 


2.36±0.08 


2.00± 0.069 


Arg(5) (deg) 


15.1 ± 3.5 


11.2 ± 3.9 


12.2 ± 3.5 


12.5 ± 3.1 


12.6 ± 3.1 


P 


-0.137±0.003 




-0.079±0.003 






O-HLS 




2.3690±0.011 








xVdof 


34.4/24 


18.1/24 


18.6/24 


39.1/25 


39.4/25 



dependence in the value of the phase. Note that, in contrast to the situation for the other 
data sets, only the HLS and GS models provide acceptable quality fits to the KLOE data. 
That the phase error is so much greater for the KLOE(02) fits than for any of the other 
data sets is likely a reflection of the small number of data points (five) in the interference 
region (770-800 MeV) compared to CMD-2(94), CMD-2(98), and SND which reported, 
respectively, 10, 12, and 13 data points in that range and so provide a significantly better 
characterization of the interference shoulder. 

Looking across models within each data set, we find that the HLS, GS, GP/CEN"*", 
and GP/CEN"'"'" models favour a slightly larger p mass than does the KS model. For 
comparison. Belle has recently reported the analogous values for the charged p parame- 
ters, mp± = 774.6 ± 0.5 MeV and Tp± = 148.1 ± 1.7 MeV, obtained from vr^^vrV^ 
using the GS model for the p, p' and p" Breit-Wigner shapes [33| . 

With three exceptions, the quahty of the fits, as measured by the x^/dof, is reasonable. 
The fits to the KLOE(02) data using the KS, GP/CEN+ and GP/CEN++ models yield 
unacceptably large x^/dof, and hence are not included when arriving at a final combined 
assessment of our results below. The GS and HLS models provide the best quality fits 
in all cases, and the only acceptable quality fits to the KLOE data. The fit residuals 
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TABLE III: Results of fits to the SND data. 



Parameter 


KS 


HLS 


GS 


GP/CEN+ 


GP/CEN++ 


nip (MeV) 


771.60±0.4d 


773.42ib0.45 


774.26±0.47 


776.37ib0.44 


776.37ib0.44 


Tp (MeV) 


142.77±0.78 


147.23±0.89 


147.34± 0.84 


(147.34) 


(147.31) 


5Tp (MeV) 








-0.32ib0.42 


-0.35ib 0.42 


\S\ (10-3) 


1.96±0.07 


2.04±0.07 


2.02±0.07 


2.47ib0.08 


2.09ib 0.07 


Arg(5) (deg) 


16.2 ± 1.8 


15.1 ± 1.7 


15.4 ± 1.7 


15.8 lb 1.6 


15.9 lb 1.6 




-0.144±0.003 




-0.081±0.003 










2.3770±0.0098 










34.5/40 


23.6/40 


25.1/40 


62.6/41 


62.9/41 



TABLE IV: Results of fits to the KLOE(02) data. 



Parameter 


KS 


HLS 


GS 


GP/CEN+ 


GP/CEN++ 


rup (MeV) 


771.21ib0.18 


773.05ib0.17 


773.90ib0.18 


776.91ib0.16 


776.91ib0.16 


Tp (MeV) 


141.88ib0.36 


147.29ib0.41 


147.04ib 0.39 


(146.50) 


(146.46) 


5Tp (MeV) 








-2.29ib0.21 


-2.33ib 0.21 


|<5| (10-3) 


1.26ib0.08 


1.49ib0.08 


1.48ib0.08 


2.32ib0.10 


1.98± 0.08 


Arg((5) (deg) 


19.4 lb 7.4 


6.0 lb 6.3 


9.4 lb 6.3 


23.2 lb 4.5 


23.5 lb 4.5 




-0.145ib0.0011 




-0.085±0.0011 






auLS 




2.3880ib0.0034 








xVdof 


293.2/55 


63.6/55 


72.4/55 


315/56 


318/56 



(not shown here) reveal that the GP/CEN"*" and GP/CEN"'"+ models tend to lie below 
the data above ^/s = rUp, which suggests that, with the high quality of the current data 
sets, forms of the resonance chiral effective theory model which include explicit higher 
resonance (in particular p') contributions will be needed if one wants to further improve 
the fit quality. 



IV. IMPLICATIONS FOR a, 



As is well known, the leading order hadronic vacuum polarization contribution to 
can be obtained from the experimental e~^e~ — )■ hadrons cross section using the dispersive 
representation 



14 



^LO 
VJ had 



a 



EM 



(0) 



37r2 



4ml ^ 



(16) 



where K{s) has the form given in [IJ] and R{s) is the ratio of the "bare" e~^e~ — )■ hadrons 
cross-section to that for e~^e~ — ?• pC^ pT . The integrand in Eq. [161 terms of F.,^{s), has 
the form R{s) = (1 — 4m^/s)3/^|F7r(s)p/4. Fitting the data using any of the models for 
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TABLE V: [(^a/,]^fd;m« ^ 10^° ^^r the models discussed in the text and the CMD-2, SND, and 
KLOE e^e~ — )• tt^tt" cross-sections. 



Experiment 


KS 


HLS 


GS 


GP/CEN+ 


GP/CEN++ 


CMD-2(94) 
CMD-2(98) 

SND 
KLOE(02) 


3.8 ±0.6 
4.0 ±0.6 
4.2 ±0.4 
(2.2 ±0.6) 


4.0 ±0.6 
4.6 ±0.6 
4.3 ±0.4 
4.2 ±0.7 


2.0 ±0.5 
2.5 ±0.5 
2.2 ±0.3 
2.2 ±0.6 


2.0 ±0.5 
2.2 ±0.4 
1.9 ±0.3 
(0.5 ±0.8) 


1.8 ±0.4 
2.1 ±0.4 
1.7 ±0.3 
(0.3 ±0.8) 



F^(s) discussed above, in which the IB contribution is that associated with the u, thus 
allows a straightforward, albeit model-dependent, separation of the p-u IB and nominally 
IC components of the experimental cross-sections. Writing Ft,{s) = Ft,^nic{s) + F.,,^rnix{s) , 
with Fj^^rnixis) and -F7r,Ar/c(s) the p-u IB and nominally IC (NIC) {p+p'+- ■ ■) contributions 
to Ft,{s), and defining (5|F^(s)p = |F^(s)p — I-Ftt.at/cI-s)!^, we obtain the desired p — u 
mixing plus direct — ?■ tttt coupling contribution to [^ct/J/j^d' 

[^«.]L%.. = ds^il - 4ml/s)y^6\FAs)\\ (17) 

This represents one component of the full IB correction which must be applied to the 
experimental r — )■ ttttz/^ data, the other components being those associated with (i) the 
impact of the vr^ — 7r° mass difference on the phase space factor, (ii) possible differences 
in the p° and p"^ masses and widths, (iii) possible differences in the products of decay 
constants and vrvr couplings of the various charged and neutral p resonances, and (iv) 



long-distance electromagnetic corrections |36H40| . Many, though not all, of recent r- 
based analyses of [a^ J^^ |20l. 23-27, 34j_35[ employ a value of [Sa^]had„iix obtained using 



the GP/CEN model [36|, l37[. In Ref. [4l| it was pointed out that the generic structure of 
the p — uj interference contribution to Ft^{s) and the monotonically decreasing nature of 
the integration kernel, K{s), combine to introduce strong fit-parameter-sensitive cancel- 
lations in the integral in Eq. [T71 and hence significant model dependence into [^ci^^had rnix 
(see below for more on this point). The variation in the values of [5a/i]^^.mia; across the 
various models was found to be greater than the experimental uncertainty produced by 
any single model, and to necessitate a re-evaluation of the overall uncertainty on this 
quantity. 

Updating the analysis of Ref. jiH to incorporate the data improvements described in 
Section Uni we obtain the values of [^cii^hadmix shown in Table IVl for the models described 
above (with s^ax = 2.25 GeV^). The results for the KLOE(02) data ± (KS, GP/CEN+, 
GP/CEN"'"+) model cases are shown in brackets to highlight the fact that the underlying 
fits are of poor quality and that these values are not used in fixing any of the central 
values, or error ranges, quoted below. 

We see that, despite the higher statistics of the CMD-2 (98) data, the uncertainty in 
i^^f^had-mix has uot been reduced. This is because of the increased uncertainty in the 
phase of the IB term resulting from the poorer beam energy calibration. The corrections 
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to the original SND data (the SND results quoted in Ref. [4l|] were based on the un- 
corrected results) turn out to have little impact on [Sa^]j^^^.^^^. The introduction of the 
new KLOE(02) data and covariance information allows us to significantly improve on our 
earlier treatment of the now-superceded KLOE(Ol) results. We note, in particular, that 
the KLOE(02)-based values for the HLS and GS models are now consistent with those 
obtained from the other two experiments. 

In arriving at a final assessment of our results for 0,11] had- mix ^ have adopted the 
view that, since all the models considered have a reasonable basis in phenomenology, all 
results corresponding to a given data set and given model which produce an acceptable 
quality fit are to be included in the assessment. Our final results (both here, and for the 
quantities to be discussed in the next section) are thus obtained by first performing a 
weighted average over all experiments for each separate model, and then taking the av- 
erage (half the difference) of the maximum and minimum values allowed by the resulting 
error intervals for the different models to define our central values (model-dependence- 
induced uncertainties). The result of this prescription for [^a^il^^.^j^, is 

[S^^^Z■,m^x = (3-1 ± l-5we/ ± 0.3,,*,) X 10-1°, (18) 



whose central value is compatible with that reported in Eq. 15 of Ref. [41 



Two features of the analysis worth emphasizing are the strong sensitivity of [ScLfi]had-mix 
to the phase, 0, of the IB contribution and the resulting sensitivity to the model employed 
for the broad p contribution. Both can be understood by considering the generic form of 
the flavor '38' component of the EM cross-section. 

Writing the isospin-breaking amplitude in the "p — u mixing" form, generically oc 
Bp{s)6Puj{s), with 6 = l^le**^ and -Bp(s) the model-dependent form of the p(770) Breit- 
Wigner, the IB, flavor '38' component of the EM cross-section is oc 



\B,is)\' 



[{ml - s)2 + mlTl] 



(cos(0)(m^ - s) - m^T^ sin(0)). (19) 



The coefficient multiplying cos(0) is antisymmetric in s about m^, which means that the 
contribution of this term to [S(ifj]had;mix vanishes in the limit that the s-dependence of 
Bp{s)K{s)/s is neglected. The result is a relative enhancement of the contribution of 
the sin(0) term, and so a significant reduction of [So,p]j^^^.^^^ relative to what would be 
obtained for </> = 0, for even modest values of (p. Thus, relaxing the = assumption of 
the GP / GEN model, the non-zero phase preferred by the data for the GP / GEN"*" and 
GP / GEN"*""*" models leads to the reduction seen in Table |V] as compared to the GP / GEN 
result [ScLplhad-mix = (3.5 ± 0.6) x 10"^ which was obtained from the uncorrected GMD- 



2(94) data, and fixed to [3J, l37|, |57 



Regarding the sensitivity to the model chosen for the form of the p contribution, the 
~ 1.5 X 10^1° model-dependence uncertainty noted above is significantly greater than 
the ~ 0.5 X 10^1° uncertainty generated by experimental errors for any given model. 
For a given model, however, the different data sets produce results which are in good 
agreement. 
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Before continuing, it is worth commenting on one feature of the relation of our results 
to those of Ref. [2l|. The latter reference, employing the KS' model, finds results for 
\^'^tj\had-mix Compatible with those obtained using the GS model, in contrast to our results 
above. The reason for this difference is that we employ the KS model, whose s-dependence 
in the region of the interference shoulder is slightly different from that of the KS' model. 
We have verified that the KS' model indeed produces results compatible with those of the 



GS model, as found in Ref. [2l|, the values corresponding to the CMD2(94), CMD2(98) 
and SND(06) data sets being 2.1(5) x 10"^°, 2.2(6) x 10"^° and 2.3(3) x 10"^°, respectively. 
The difference between the KS and KS' model results serves as a further illustration of our 
point about the model sensitivity of the determination of [^c^'^^^d■m^x■ This sensitivity is a 
consequence of the strong cancellation in the integral for [^o,fj\^d-mix-> ^ cancellation which 
enhances the impact of small differences in the s-dependence of the Fj^^mixis) among the 
different models. 

With the KS model form, the IB parameter 6 also occurs in the nominally IC p con- 
tribution to the amplitude. One might think it appropriate to also remove the effect of 
this IB contribution in comparing the r and electroproduction results. Such a procedure, 
however, would involve double counting when used in combination with a separate as- 
sessment (such as that discussed in Refs. [3^, [13]) of the IB contribution to the broad p 
component of the ratio of charged and neutral current versions of F^. The reason is that 
any 5-dependence in F^^at/c* represents only a partial assessment of IB in the coefficients 
1/[(1 + 6){1 + /3 + 7)], '/3/(l + /3 + 7) and 7/(1 + /3 + 7) of the p, p', p" contributions 
to -F^(s), coefficients which are proportional to the products fpgp-jm, fp'Op'-Kn, ■ ■ ■ of the 
various p resonance decay constants and vrvr couplings. An assessment of IB contributions 
to these products is, however, already provided by the treatment of Refs. [36|, 



While the p—u interference shoulder, because of its narrow structure, is unambiguously 
identifiable as an IB effect associated with the intermediate u state, there is no signal, 
internal to the electroproduction data, that allows one to identify IB contributions to the 
broad p, p', p" parts of the cross-section. Thus while, for the KS model, the presence of 
the IB parameter 6 in the broad p contribution produces an explicit IB contribution to 
fpOpirn, additional implicit IB contributions are, in general, also unavoidably present in 
the coefficients + /3 + 7), f3/{l + 13 + j) and ■y/{l + (3 + 7). Indeed, if one simply 
reparametrizes the KS model, writing 

^(KS)(^) ^ I^ Pp(g) (1 + 6P^{s)) + P'PAs) + yPp"(^) ^ ^ (20) 

explicit IB dependence, which in the original parametrization was present only in the p 
and p — u "mixing" contributions, is now present in all of the p, p', p" and p — u "mixing" 
contributions. The presence or absence of such explicit IB in the nominally IC p, p', p" 
contributions to Ft^{s) is thus parametrization-dependent, from which it follows that the 
integrated contribution associated with such explicit IB contributions for any particular 
parametrization has no well-defined physical meaning. 

While framed in the context of the KS model, the above discussion, of course, applies 
more generally; implicit IB contributions, which cannot be identified experimentally, will 
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also be present in the NIC part of the fitted version of F^^ obtained using any model. For- 
tunately, there is no need to tackle the difficult task of trying to identify such additional 



implicit IB contributions; Refs. [36|, |37| provide a theoretically well-founded framework 
for assessing the full set of IB contributions to the NIC part of F^, one that is expected to 
produce reliable results in the region below s ~ 1 GeV^ which dominates the SM contri- 
butions to [a^]^^- With sufficiently good r decay and electroproduction data one could, 
of course, instead, fit any particular model separately to both data sets and identify 
small IB effects in the fitted NIC contributions that result. Current data is insufficiently 
precise to make such an experimental alternative feasible at present. 

We note finally that, independent of the above discussion, one could create a third 
variant of "the" KS model by replacing the factor s/m^ in the mixing term of the KS' 
model with a function f{s) which varies smoothly from at s = to 1 in the region 
around the u. Such a variant would share with the KS' version of the model the absence of 
explicit IB in -F^,Ar/c(s), but produce a numerical result for the p — u "mixing" correction 
matching that of the KS version. 

As previously mentioned, the p — u mixing contribution is a part of the ratio of EM 
to weak form factors, which is itself one of three components of the overall IB correction. 



Rib, to the r-based determination of [Sa^]f^^^ considered in Ref. [37|. The uncertainty in 



the contribution of the ratio of form factors also happens to be the dominant uncertainty 
in the total correction [Sa^]^^^ associated with Rib- The results presented in Eg. [T8] show 
that the p — u contribution to the form factor ratio has a significant model-dependence, 
necessitating an increase of ~ 1.5 x 10~^° in the overall uncertainty assigned to [^a^t]^^- 
We argue that such model dependence should be added linearly to the remaining uncer- 
tainty to refiect the range that would be obtained by using different phenomenologically 
acceptable model fits. The overall uncertainty in the ratio of form factors thus shifts 
up to ±3.9 X 10~^° if one uses the figures of Ref. js^] for the other IB corrections, and 
up to ±3.1 X 10^^'^ if one uses instead the more recent results of Ref. 2l| for these cor- 



rections. Modifying the r-based assessment of Ref. [21| to account for both the slight 



shift represented by our version of the central value and our increased mo del- dependence 
of [So,fi\had-m.ix^ obtaiu the following modified version of the difference between the 
experimental and r-based SM assessment of 



a^^P -al = (14.5 ± 9.7) x 10"^° (21) 

where, as above, we have added the contribution associated with the increased model- 
dependence of [^O'fjl^dmix li^isarly to the other errors. 



V. EXTRACTION OF Hp,^ AND THE DIRECT w ^ tttt COUPLING 

A. Analysis 

In general it is unavoidable that, for a given choice of p and u interpolating fields, 
the p — u interference signal in F-,^{s) will be the result of a combination of true p — u 
mixing and a direct (non-mixing-induced) IB a; — )■ tttt contribution. It is of interest to 
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separate these contributions in order to gauge the strength of the intrinsic wtttt couphng 
and to estimate the off-diagonal p — uj element, np(^(g^), of the vector meson self-energy 
matrix. Such a separation is relevant for meson-exchange models of IB in the NN inter- 
action, which models require information about np^(g^) at < 0. Often an "effective" 
^pujil'^), extracted from electroproduction data under the assumption that direct u rni 
contributions can be neglected, has been used for this purpose (see, e.g., Ref. [58|). In 
general, such an estimate will be contaminated by direct cu — >■ tttt contributions, which 
are unavoidably present in the electroproduction cross-section in the interference region. 

The separation of mixing and direct a; —t- tttt contributions depends on the model 
used to represent the broad p contribution to -F^(s). The feasibility of such a separation 



was discussed in detail in Ref. [12(], though the lower quality of the data then available 
0] allowed only weak constraints to be placed on IIp^ and the strength of the direct IB 
a; — )■ TTTT coupling, (/^/^Tr. Below, we revisit this analysis using the recent data described 
in section UTTl and the various models of Ft^{s) described in section ITTl 

In the interference region, the time-like EM pion form factor is given, in terms of the 
physical vector meson degrees of freedom, by 



FJs) 



+ background 



(22) 



where Dyv{s) are the scalar propagators of the physical vector mesons, gp.„T, and g^T.^, 
are the couplings of the physical p and uj to the two-pion final state, and /^^ and /^^ are 
the electromagnetic p and uj couplings. The background term includes all non-resonant 
contributions to the form factor. 

The physical p and uj fields referred to in Eq. ( !22l) are admixtures of the isospin-pure 
1 = 1 and / = fields and can be written as 



P = P^ - (^1^^, 



UJ 



(23) 



with ei and 62 as given in Ref. [12] to first order in isospin breaking. The physical 
UJ — )■ TTTT coupling, (?aj7r7r, thus coutaius both a mixing contribution and a direct IB coupling 
contribution. Explicitly (l2| . 



9u 



9^1^^'^ - z)- izTg 



(24) 



where 



z = 



rhpTp 



— t 



rhpTp 



(25) 
(26) 



with Ilp^ the real part of flp^j and my the real part of the complex pole position of the 



vector meson V. The imaginary part of IIp^ is given by 



-GrhpTp with G 



In the limit that and m^— are negligible, z = 1, and the first term in Eq. [24lvanishes 
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identically, making a determination of g^i t^t^ impossible in this limit. Fortunately, z turns 
out to deviate sufficiently from 1 to allow a reasonable determination of gc^ij^T^- 
With the above definitions, Eq. [22] becomes 



Pp{s) + |rex|e^'^=+=- ((1 - z)G - izT)P^{s) + background (27) 



where = rex = |rex|e*''e+e- and, experimentally [59[, 



rh?S{p — e+e^ 



.1/2 



0.296 ±0.005 



(28) 



using PDG07 [60] values for the widths and masses. The full propagators appearing in 
Eq. [22] have been replaced in Eq. [27] with the same forms as were used in Sec. [TT] (Eq. [2] for 
the KS and HLS models, Eq. [8]for the GS model, and Eq. [TO]for the GP/CEN family), 
but with their dimensionful numerators removed (e.g., Pp{s) = Pp{s)/m?p for the KS and 
HLS models). (This preserves the usual dimensionful definition of /y^.) 

Fits to experimental data provide only the magnitude of the u contribution to Ft^{s) 
and its phase relative to the dominant p term, i.e., the modulus and phase of the coeffi- 
cient of Pui{s) inside the square brackets in Eq. [27] This coefficient is written compactly 
as Ae^'^ with the so-called Orsay phase. It follows from the above that the values of 
A and (f) obtained from a given data set depend on the precise form of the p propagator 
in the model used to represent Ft^{s). Any separation of G and T is thus unavoidably 
model-dependent. 



B. Numerical Results 

The Orsay amplitude and phase corresponding to the data sets and models described 
in sections [TT] and IIIII are shown in the first two rows of Tables IVII to IIXI The third 
and fourth rows of these tables also show the real and imaginary parts of z, which 
are determined in each case by the corresponding fit values of the p mass and width. 
The range of Orsa y p hases found here is consistent within errors with the old value of 



(j) = 116.7° ± 5.8° [59j obtained using less precise data and a model of Ft,[s) similar to 
the KS model but with a constant background term instead of the p' resonance. The 
central value of our result for the phase, obtained as described previously, is then 

0=113° ±4;;,,,,, ±2°,,, (29) 



slightly lower than the result of Ref. [59[, but consistent with it within errors. Omission 
of the KLOE data makes no significant difference to this determination. 

The Orsay phase uncertainty for any single model ranges from 2° for the high-precision 
SND data to 5° for the KLOE(02) data. As expected, the phase determination is slightly 
less precise for the CMD-2(98) data than for the lower statistics CMD-2(94) data because 
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TABLE VI: Orsay amplitude, phase, and separated mixing and direct wtttt coupling parameters 
for the CMD-2(94) data. 



Parameter 


KS 


HLS 


GS 


GP/CEN+ 


GP/CEN++ 


A 


0.0101 ±0.0008 


0.0103 ±0.0008 


0.0103 ± 0.0008 


0.0107 ± 0.0008 


0.0107 ±0.0005 


(f) (deg) 


114 ±3 


112 ±3 


111 ±3 


110 ±3 


110 ±3 


Re{z) 


1.041 ±0.003 


1.046 ± 0.002 


1.048 ± 0.002 


1.052 ±0.002 


1.053 ±0.002 


Im(2;) 


0.158 ±0.009 


0.133 ± 0.009 


0.121 ±0.009 


0.102 ±0.009 


0.103 ±0.009 


G 


0.080 ± 0.032 


0.073 ± 0.032 


0.076 ± 0.032 


0.074 ± 0.032 


0.074 ±0.031 


f 


-0.043 ± 0.005 


-0.042 ± 0.004 


-0.041 ± 0.004 


-0.041 ± 0.003 


-0.041 ± 0.003 


sin (/)e+e- 


0.120 ±0.004 


0.124 ±0.003 


0.124 ±0.003 


0.129 ±0.003 


0.129 ±0.003 



TABLE VII: Orsay amplitude, phase, and separated mixing and direct utttt coupling parameters 
for the CMD-2(98) data. 



Parameter 


KS 


HLS 


GS 


GP/CEN+ 


GP/CEN++ 


A 


0.0107 ± 0.0004 


0.0107 ±0.0004 


0.0107 ± 0.0004 


0.0111 ±0.0004 


0.0111 ±0.0002 


(j) (deg) 


114 ±4 


109 ±3 


109 ±4 


108 ±4 


108 ±3 


Re(z) 


1.041 ±0.002 


1.046 ± 0.002 


1.049 ± 0.002 


1.054 ±0.001 


1.055 ±0.001 


Im(z) 


0.153 ±0.008 


0.128 ± 0.008 


0.116 ±0.008 


0.087 ± 0.008 


0.087 ±0.008 


G 


0.083 ± 0.037 


0.044 ± 0.040 


0.055 ± 0.037 


0.058 ±0.036 


0.058 ± 0.034 


f 


-0.046 ± 0.005 


-0.039 ± 0.004 


-0.040 ± 0.004 


-0.040 ± 0.003 


-0.040 ± 0.003 


sin 


0.127 ±0.003 


0.130 ±0.002 


0.130 ±0.003 


0.135 ±0.003 


0.135 ±0.003 



of the former's larger correlated systematic errors. The precision of the phase determi- 
nation from the new data, for a given model, has been improved by a factor of ~ 1.5 — 2 
compared to that from the earlier data. 

The 30% to 50% reduction in phase uncertainty significantly improves the reliability 
of the separation of pure mixing and direct ujittt coupling contributions to Ftj{s). The 
results for G and T corresponding to the best fit values for A and are shown in rows 
5 and 6 of Tables |VT] to IIXI The one-sigma region in the G-T space, shown for_ each 
data set and model in Figures [T] |2l |3l and HI defines the uncertainties in G and T and 
is based on full propagation of the fit parameter covariance matrix, itself a reflection of 
the experimental covariances. 

As Figures [T]-|l] show, it is possible to constrain T rather well. With the prescription 
above for handling averages and uncertainties, we obtain 

f = -0.044 ± 0.006„orfe/ ± 0.002rf,ta. 

The consistency across different models for each data set is actually quite good. 

The determination of G is less precise, with the prescription above for determining 
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TABLE VIII: Orsay amplitude, phase, and separated mixing and direct utttt coupling param- 
eters for for the SND06 data. 



Parameter 


KS 


HLS 


GS 


GP/CEN+ 


GP/CEN++ 


A 


0.0110 ±0.0004 


0.0112 ±0.0004 


0.0111 ± 0.0004 


0.0116 ± 0.0004 


0.0116 ±0.0002 


4> (deg) 


116 ±2 


113 ±2 


113 ±2 


112 ±2 


111 ±2 


Re{z) 


1.037 ±0.002 


1.044 ± 0.002 


1.047 ± 0.002 


1.053 ±0.001 


1.053 ±0.001 


Im(2;) 


0.166 ±0.007 


0.136 ± 0.007 


0.124 ±0.007 


0.094 ± 0.007 


0.094 ±0.007 


G 


0.097 ±0.020 


0.087 ±0.020 


0.090 ± 0.020 


0.093 ± 0.020 


0.094 ±0.019 


f 


-0.050 ± 0.003 


-0.046 ± 0.003 


-0.045 ± 0.002 


-0.044 ± 0.002 


-0.044 ± 0.002 


sin (/)e+e- 


0.130 ±0.003 


0.133 ±0.003 


0.133 ±0.003 


0.140 ±0.003 


0.140 ±0.003 



TABLE IX: Orsay amplitude, phase, and separated mixing and direct lotttt coupling parameters 
for the KLOE(02) data. 



Parameter 


KS 


HLS 


GS 


GP/CEN+ 


GP/CEN++ 


A 


0.0071 ± 0.0003 


0.0081 ± 0.0003 


0.0081 ± 0.0003 


0.0111 ±0.0003 


0.0111 ±0.0003 


4> (deg) 


119 ±6 


104 ±5 


107 ±5 


118 ± 5 


119 ±5 


Re{z) 


1.034 ±0.001 


1.042 ± 0.001 


1.046 ± 0.001 


1.055 ±0.0006 


1.055 ± 0.0006 


Im(z) 


0.173 ±0.004 


0.141 ± 0.004 


0.129 ±0.004 


0.087 ± 0.003 


0.088 ± 0.009 


G 


0.100 ±0.041 


0.006 ± 0.042 


0.028 ± 0.042 


0.169 ±0.054 


0.173 ±0.050 


f 


-0.038 ± 0.006 


-0.027 ±0.005 


-0.030 ± 0.005 


-0.047 ±0.004 


-0.047 ±0.002 


sin 


0.080 ± 0.002 


0.100 ±0.001 


0.099 ± 0.001 


0.129 ±0.003 


0.128 ±0.003 



the average across models and mo del- dependence uncertainty yielding 

G = 0.080 ± OmQ^odei ± 0m5data. 

For comparison, the uncorrected SND analysis [5i] , using a model very similar in form to 
the KS model described in Section|Tll reported G = O.lliO.Ol, consistent with our result, 
G = 0.10±0.02, for the SND±KS case. The present analysis favors a non- vanishing G, as 
expected on general grounds, with a vanishing G lying ~ 2a away from our central value. 
An improved determination of the Orsay phase, particularly in the case of the KLOE 
data, is required to make further progress on this issue. To illustrate this point, and also 
because of the lingering KLOE shape discrepancy [IS], if one omits the KLOE data from 
the experimental averages one obtains results with a slightly reduced model-dependence 
uncertainty, 

T = -0.045 ± O.OOA^odel ± 0m2data 

and 

G = 0.084 ± 0.022^odei ± 0.016^,^,. 

The values for G and T obtained in Tables IVIHIXI imply, using the earlier definitions, 
and our prescription for combining the results corresponding to different models, the 
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FIG. 1: CMD-2(94) G and T one-sigma regions. 



result 

fipclmj) = -4950 ± 530 Wei ± '^^^data MeV^ (30) 

for the real part of np^(m^). The imaginary part of np^(mp) is, similarly, found to be 

—8600 ± 2510mode/ ± 1690rfata MeV^. The result is an estimate of the off-diagonal part of 
the physical p — u self-energy matrix at the p mass 

np^(mj) = (-4950 ± 780) + (-8600 ± 4200)i MeV^ (31) 

where we have combined the model-dependence and experimental errors linearly. Mod- 
estly reduced uncertainties are again obtained if the KLOE(02) data is excluded in form- 
ing the combined result. Explicitly: 

np^(mp) = (-4950 ± 670) + (-8890 ± 4000)i MeVl (32) 

This quantity is not to be confused with the sometimes quoted 'effective' p — oj mixing 
matrix element, which contains direct w — t- tttt contributions. 

The separation described here is possible only because the two contributions enter 
the form factor with different (but known) phases, as shown by Eq. |271 Experimental 
determination of the combined phase (the Orsay phase) then allows the magnitude of each 
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FIG. 2: CMD-2(98) G and T one-sigma regions. 



contribution to be determined separately. Any improvement in the separation of the p—u 
mixing and direct a; — )■ vrvr IB contributions thus requires corresponding improvements 
in the determination of the Orsay phase. 



VI. CONCLUSIONS 



Recent experiments provide a much-improved characterization of the prominent p — oj 
interference shoulder in the e^e~ — )■ 7r"'"7r~ cross-section, making precision analyses of 
isospin breaking in the p — u system possible. Such an analysis is inherently dependent 
on the model chosen to represent the broad p resonance contribution. By studying a 
range of models, all of which provide good fits to the data, we have shown that the 
model sensitivity of most fit parameters is generally greater than the statistical uncer- 
tainty associated with any given model. The only exception is the phase, Arg(5), of 
the parameter describing the strength of the u contribution, for which the statistical 
uncertainty tends to be greater than the model sensitivity. 

Any model sensitivity of course propagates to all quantities obtained from the fit 
parameters. The example discussed in Section IIVI is that of the p — u interference con- 
tribution to the IB correction, [(^a^]^^, required for incorporating hadronic r decay data 
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SND 



-0.025 




FIG. 3: SND G and T one-sigma regions. 



into the determination of the SM expectation for [a^]j^^^. Using our result, Eq. fTSlthe 
full IB correction obtained in Ref. js^l using the method of Ref. [s^] (see alsojSTjl ) . 
Wh^d = (-13.8 ±2.4) X 10-1° (or [Sa^]^^^^ = (-16.6 ± 1.6) x 10-^° based on the recent 
analysis of Ref. |2l|), receives only a small shift of —0.4 x 10^^° (+0.7 x 10"^°) in its 
central value. The uncertainty in the size of the correction, however, increases by a much 
more significant 1.5 x 10"^°. In arriving at the combined error, we have added the com- 
ponent associated with the mo del- dependence of the p — u mixing contribution linearly 
to the other errors, for the reason discussed above. While the increased error does not 
serve to resolve the ~ (17 ± 10) x 10"^° [Tsl, Hoj (~ (11 ± 5) x 10"^° (21)) r vs. EM 



discrepancy which existed prior to the preliminary BaBar ISR results, it is nonetheless 
important to take into account in assessing the uncertainty in the IB corrections which 
must be apphed if one wishes to use r decay data in evaluating the SM contribution to 



\a 



V-ihad- 



The precision of the recent electroproduction data also allows new improved results for 
the (model-dependent) separation of contributions to Ft^{s) from p — u mixing and direct 
(non-mixing) IB w — > tttt transitions. We estimate the ratio of the two-pion couplings of 
the isospin-pure u and p to be G = 0.080 ± 0.041 where the uncertainty is conservatively 
chosen to encompass the experimental error intervals of all of the models considered. 
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FIG. 4: KLOE G and T one-sigma regions. 



Further improvements in the determination of G may be possible in future. As can be 
seen in Figures [T]-|ll the off-diagonal part of the vector meson self-energy matrix evaluated 
at the p mass is very well constrained by the analysis. The central value we obtain for the 
real part is given in Eq. [301 where the uncertainty again conservatively encompasses the 
full ranges corresponding to all models and all data sets. Our central value of Ilpi^{rn?p) 
which is free of contamination from direct uj ^ tttt coupling contributions, and thus differs 
from the usually quoted 'effective' p — u mixing matrix element obtained by ignoring 
such a direct coupling, should be used in place of the 'effective' matrix element results 
in meson-exchange models of IB NN scattering. 

Further improvement of the separation of direct and mixing contributions requires a 
reduction in the uncertainty of the Orsay phase, which in turn requires a further increase 
in the precision of the e^e^ — > tt+tt^ data in the p — u interference region. This requires 
not only increased statistics but also very precise beam energy calibration such as that 
expected from VEPP-2000. 

NOTE ADDED: Subsequent to the posting of this paper, the BaBar collaboration posted 
its results for the e+e~ — t- tt+tt" cross-sections, obtained using the ISR method Wi l. We 
will report the results of the analysis discussed here, applied to that data set, as soon as 
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the data becomes publicly available. 
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